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The so-called piv puzzle of J/ip and tp{2S) decays is examined using the experimental data avail- 
able to date. Two different approaches were taken to estimate the ratio of J/tp and ip{2S) hadronic 
^~> , decay rates. While one of the estimates could not yield the exact ratio of ip{2S) to J/tp inclusive 

("^ ■ hadronic decay rates, the other, based on a computation of the inclusive ggg decay rate for ip{2S) 

f^ ' {J I'i') by subtracting other decay rates from the total decay rate, differs by two standard deviations 

CN| . from the naive prediction of perturbative QCD, even though its central value is nearly twice as large 

s^ ' as what was naively expected. A comparison between this ratio, upon making corrections for specific 

Q I exclusive two-body decay modes, and the corresponding experimental data confirms the puzzles in 

J/i/> and ip{2S) decays. We find from our analysis that the exclusively reconstructed hadronic decays 
of the 'ip{2S) account for only a small fraction of its total decays, and a ratio exceeding the above 
C^^ , estimate should be expected to occur for a considerable number of the remaining decay channels. 

Cn ■ We also show that the recent new results from the BES experiment provide crucial tests of various 

theoretical models proposed to explain the puzzle. 

m ■ 

>\ 
\o . 

O; 

>! One of the outstanding problems in heavy quarkonium physics is the strong suppression of the ip{2S) decays to vector 

,__! plus pseudoscalar-meson (VP) final states, pn and K^K + c.c, which is referred to as the pir puzzle [|l[. Following 

^\ the first observation of this anomaly g], meagre experimental progress was made over the years, and theoretical 
0^ ' analysis based on limited data often lead to unsatisfactory, sometimes premature, inferences. The situation has been 
r^ changed dramatically in the last few years. A wealth of interesting new information, which extended the puzzle 

^ 1^ considerably, has emerged from intense studies of '0(25') hadronic decays at the BES experiment, using a large sample 
JL 1 of 3.79 million '0(25') decays 0. It is hoped for that new concerted efforts on both the theoretical and experimental 
^ I side would eventually lead to a solution of this long-standing conundrum. 

r^ . In this paper we seek to examine the pir puzzle based purely on existing experimental data. We begin with an 

Kjl ■ analysis for estimating the ratio of hadronic decay rates of J /ip and '0(25), which we shall denote by Q, by using the 

data compiled by the Particle Data Group B, in an attempt to avoid as many theoretical ambiguities as possible 

in the analysis. Two different approaches to this estimate are performed. First we compare the results between 

j^ , themselves, and then the naive prediction of perturbative QCD (PQCD) is used. Subsequently, possible corrections 

■ " ' to Q are discussed, as they associate with specific exclusive decay modes, and the corrected values of Q are used as 

standards to compare with the corresponding experimental data. Comments on the issue of J/'0 and '0(25) decays to 

multi-hadron final states and on the potential similarity of the rjc-rjci'^S) decays to the pn puzzle are profusely added. 

Finally, various theoretical models are discussed in order to offer a greater understanding of the pir puzzle, in light of 

recent BES results. 

Conventionally, measured ratios of ipi'^^) to J/0 branching fractions for specific exclusive hadronic decays are 
compared with the naive prediction of PQCD, the so-called "15% rule". In the framework of PQCD [g, one expects 
J/ip{'ip{2S)) to decay to hadrons via three gluons, or a single direct photon. In either case, the partial width of the 
decay is proportional to |\['(0)|^, where 1^1/(0)1 is the wave function at the origin in the nonrelativistic quark model of 
the cc quark state. Thus one finds that 

_ BW2S) -. ggg) _ a^(V^(25)) i3(0(25) ^ e+g-) _ ,, , , ^ , ,,o, ... 

^' = B{J/i.^ggg) " a^^J/^) S(J/V. - e+e") " ^^^'^ ^ '■'^^"' ^'^ 

where the new world averages of the leptonic branching fractions are used (see Table 1). This is assuming that the 
strong coupling constants are equivalent, i.e., as{ip{2S))=as{J/'4')- Taking the running constant as into account [||, 
the ratio Qh becomes 

Q,, = (12.5 ±1.9)%. (2) 

The MARKII experiment first compared the theoretical prediction of this value ((12.2 ± 2.4)% then used) with 
measurements for a number of exclusive hadronic decays of the J/tp and the ip{2S), thus revealing the pn puzzle [^. 



However, this naive prediction suffers several apparent approximations. Higher order corrections, which may not 
even be small, are not included in this calculation. For example, a first order correction to the branching fraction 
of J/ip -^ e+e" could be 50% of the lowest term if one were to use as{mj/^) '--^ 0.2 H. The relativistic effect is 
also ignored. Since the mass difference between J/tp and V'(2S') is around 20% and < tr /c^ >'^ 0.24 for J/f/", this 
correction may be at the same level as the lowest order |0] . The inclusion of the finite size of the decay vertex will 
significantly reduce the ggg decay width of J/^ j^]. Moreover, the effect of non-perturbative dynamics is neglected, 
the size of which is hard to estimate. Therefore, people may question the validity of the "15% rule" as a serious 
benchmark for comparing experimental data. 

We present here two approaches to estimate the ratio Qh, using the data as displayed in Table f-lH. The first 
approach is based on an assumption that the decays of the J/t/i and tp{2S) in the lowest order of QCD are classified 
into hadronic decays (ggg), electromagnetic decays (7*), radiative decays into light hadrons (735), and decays to 
lower mass charmonium states(ccX) HJl^. Thus, using the relation B{ggg) + B{'^gg) + -8(7*) + B{ccX) = 1, one 
can derive B(ggg)+B(7g5) by subtracting B(7*) and B(ccX) from unity. 

The electromagnetic decay channels of the J/ip produce hadrons, e+e~ and /^^/i^ as final states. Besides these 
channels, the electromagnetic decays of ip{1S) also include the t'^t~ as a final state. The experimental data is summa- 
rized in Table 1, where the branching fraction of ipi'^S) -^ t~^t~ is a recent measurement of the BES experiment [ [fl| , 
whereas the other data was taken from the Particle Data Group [E 12 1. The total contributions to the electromagnetic 



decays of J/^^: and ip{2S) are then given as B{J/iJj -> 7*) = (28.81 ± 2.00)% and B{i'{2S) -^ 7*) = (5.08 ± 0.55)%, 
respectively. 

As regards to the decay into lower mass charmonium states, the J/ip has only one radiative decay channel into rjc, 
whereas the ip{2S) can decay into a number of other final states, i.e. tt+tt" J/V', tt'^tt^ J/?/;, 7r° J/-0, rjj/ij:, 7XcO, 7Xcij 
7Xc2, 7?7c, lVc{2S) and l^Pi + X. The decay rates of the last two channels are faint and thus are neglected in our 
calculation. The experimental data summarized in Table 2 are all taken from PDG M. Using this data we calculated 
the total contribution to ccX: B{J/%1) -^ ccX) = (1.3 ± 0.4)% and S(V'(25) -^ ccX) = (78.1 ± 3.9)%, respectively. 

By deducting the contributions B(7*) and B(ccX), we find that B{J/ij -> ggg) + B^J/il) -^ jgg) = (69.9 ± 2.0)% 
and B{ip{2S) -^ ggg) + B{iJj{2S) -> -fgg) = (16.8 ± 3.9)%. Therefore the ratio of branching fractions of ■ip{2S) to J/^p 
decays into hadrons is given by 

n Bj^^'^S) ^ ggg) + B{i^{2S) ^ 75.9) .0. n + ^ fiW Cx^ 

^^ - BiJ/^^ggg)^BiJ/^^,gg) = ^^^'^ ± '■'^''- ^'^ 

The relation between the decay rates of ggg and 7317 is readily calculated in PQCD to the first order as 0] 

T^J/tp ^ ggg) 5 a^mc) tt 

Using as{mc) = 0.28, one can estimate T{J /ip — *■ ^gg)/T{J/ip -^ ggg) ~ 0.062. A similar relation can be deduced for 
the 4'{2S) decays. Thus one expects that "24% ratio" stands well for either ggg mode or ■^gg mode. 

The other approach is to use the data on branching fractions for hadronic decays in final states containing pions, 
kaons, and protons that have already been measured for both J/ip and ■0(25'). They are 3(7r"'"7r~)7r'^, 2(7r"'"7r~)7r'^, 
7r+7r-7r°, tt+tt-K+K-, tt+tt-PP, PP, PPtt" and K+K'. Using the PDG data compiled in Table 3, we have 

J2 B,{J/^ -^ fi) = (9.43 ± 0.72)% and Y. Bi{i}{2S) ^ fi) = (0.941 ± 0.185)%. It follows that 

8 8 

O2 = ^B.(V(25) ^ /O/^S.(J/0 ^ f^) = (10.0 ±2.1)%. (5) 

i=l i=l 

We note that the results obtained by the two approaches vary considerably. However, a comparison of the values 
of the total branching fraction for tp{2S) — > ggg computed by the two approaches indicates that only a small fraction 
{^ 6%) of the exclusive hadronic decays oi'p{2S) have been reconstructed experimentally. It is thus obvious that Q2 is 
not the exact ratio of ip{2S) to J/tp inclusive hadronic decay rates, but represents on average the ratio of the exclusive 
decay channels, as measured to date. We therefore do not consider Q2 any further. Nevertheless, the question persists 
as to where the remaining hadronic ip{2S) decay modes are and how the corresponding pattern of decays for ■0(25') 
to J /ip behaves. It would be an intriguing experiment task to search for those remaining channels that are in such 
final states as those with higher multiplicities, or those with multi-neutral particles, or even for remaining channels 
in non-gg states. In comparison with the naive PQCD expectations, the central value of Qi is about a factor of 
two higher than that of Qth, as stated in Eq. (2). However, the difference lies within the 2a error of Qi and is 



only marginally significant. The substantial error of Qi is essentially due to the propagation of errors during the 
subtraction of the decay rate B{%Ij{2S) -^ ccX) from the total rate, although the total error of B{i/j{2S) -^ ccX) itself 
merely amounts to about 5%. Taking into account the apparent approximations to the naive expectations of PQCD 
as well, it seems to us premature to regard this as a remarkable discrepancy that deserves serious considerations. 

To use the ratio Qi for comparison with the experimental results of ip{2S) and J/ijj decays, it should be noted 
that the estimate of Qi, like the prediction of Qth, is made for the total width for ggg decay, not for the partial 
widths of exclusive final states. Consequently, a number of corrections may be associated with specific exclusive decay 
modes: (I) It is shown that the J/?/; and V'(2S') decays to w/2, pa2, K*^K2 + c.c. and (/)/2(1525) are hadron helicity 
conservation (HHC) allowed ||l^, while that to pi: and K*K are HHC forbidden jlj]. The general validity of the 
HHC at the charmonium mass scale is still an open question. It is suggested that a critical test would be to measure 
the angular distributions of exclusive final states [Q . Existing measurements on angular distributions for J/^l^ decays 
into pp, AA, YPYi and K'gK^ |l^ are consistent with the HHC predictions (baryon pairs within 1 — 1.5 standard 
deviations); however, no data is available for 'ip{2S) decays so far. Exclusive processes that violate helicity conservation 
are suppressed by powers of 771^/s in QCD [Q. This would contribute a suppression factor M^,^/M3,(2S) ~ ^■'^^ ^^ 
the ratio of the 'ip{2S) to J/tp decay rates for final states pn and K*K. On the other hand, the 777 and 777' modes 



are allowed by the helicity selection rule, since helicity conservation applies only to the hadrons |13|. (II) Exclusive 



reactions which involve hadrons with quarks or gluons in higher orbital angular momentum states are suppressed 
by powers of l/.s = E^^^ Q. This contributes a suppression factor of J^'j/^/M'li2S) ^^ ^^^ ratio of i/'(2S') vs. J/%j) 

branching fractions for decays into such final states as Ci^/2, pa2, K*^K2 + c-c. and 0/2(1525), where a meson in 
P-wave state is included. (Ill) Ref. JMJ reported another suppression arising from the asymptotic form factor which 
would be M^,>/M^,2g) — 0.25 for decays to pp channel. Contrary to these calculations, Ref. n& evaluated the 
three-gluon contribution with the c-quark mass instead of the charmonium mass. As a consequence, the ratio of the 
J/ip and ip{2S) decay widths is not scaled to the 8th power of the ratio of their masses in our calculations. 

Table 4 lists corrections to the ratio Qi for several exclusive hadronic decay channels. Experimental data from 
PDG Q are also included for comparison. As seen from the table, the predicted corrected value of the ratio Qi 
for bin is consistent with the experimental data, whereas the experimental ratio of Ki{1270)K^ is enhanced as 
compared with the predicted value. The deviations of the measured ratios for decays into VP, VT and other final 
states (even pp) from the corrected ratios demonstrate suppressions, in this case. Note that the combination of all 
the above correction results lead to a substantial reduction of the ratio of ip{2S) to J/tp decay rates, well below 24% 
for many of the exclusive hadronic decay channels (this is also compatible with the observation mentioned above that 
the value of Q2 is much lower than that of Qi). One should therefore conclude that a considerable number of other 
decay channels ought to have an enhancement with a ratio above 24%, in order to make up for all these suppressed 
channels. It is puzzling that so far there has only been one channel for the ip{2S) decays observed, the K^{1270)K^ 
channel, which is enhanced relative to the J/tp |Q. Further systematic study of ■0(25') decays are anxiously awaited. 

As is seen from the above analysis, we have restricted our comparison only to decays to two-body final states. What 
if one makes a comparison for decays leading to three or more hadrons? MARKII did make such a comparison in their 
original work |2] and claimed that the ratio Q for decay modes such as ppn'^ , ppn'^'K^ , K^K^tt^tt^ , 2 (tt^^ 77^)77", and 
3(7r"'"7r^)7r*', is consistent with the naive theoretical expectations. However, it should be pointed out that most of these 
multi-hadron final states in fact include sums of several two-body intermediate states. One thus observes a mixed 
effect which may not deviate noticeably from the expected value of Q, even if a few of the two-body intermediate states 
are severely suppressed. For example, the decay J/tp -^ 2(7r"'"7r~)7r'' proceeds predominantly through intermediate 
states 6i7r, 0^/2, and a2p. The observed Q for this decay, as reported by MARKII [^, was (9.5 ± 2.7)%, does deviate, 
though not quite significantly, from the "15% rule". This results from the fact that two of these three two-body 
intermediate states are found to be anomalously suppressed p7[ |. Therefore, one must be always cautious about 
drawing conclusions from the comparison of decays of multi-hadron final states. 

An experimental situation similar to the pn puzzle occurs in the decays of the rjc in three vector meson (VV) cases, 
such as pp, K*K and (pep, and in pp. These decays are all first-order forbidden by HHC in PQCD H,EO|; however, 
they are actually observed to occur with relatively large branching fractions |J]. It is thus interesting to look for 
the analogous decays of the r]c{2S) and compare the ratio of the rjc{2S) to rjc branching fractions with the relation 
B{i]c{2S) — > ft,) ~ B{ric -^ h) predicted by Chao et al |£l|]. In testing for helicity conservation, these decay modes for 
77c and its spin-singlet partner ric{2S) play the same role as the decay modes pn and K*K do in the case of J/ip and 
its spin-triplet partner ip{2S). The search for the ric{2S) is thus important not only because the ric(2S) is one of the 
two remaining states of the charmonium family awaiting confirmation (or discovery) but also because the study of its 
hadronic decays could shed light on the puzzle of J/tp and ip{2S) decays. 



We now move on to discuss various theoretical models made to explain the pi: puzzle as it is presently formulated. 



Instead of a critical examination of many theoretical arguments (which one can find in the literature p],E2-24|), we will 
concentrate exclusively on comparing the experimental results, mostly from the BES experiment, to the predictions 
of these models, in an attempt to differentiate between them. 

The first explanation for the MARKII observation, as proposed by Hou and Soni ||2^ and generalized later by 
Brodsky et al. ^^, is the postulate that the decay J/ip, which violates the helicity selection rule of PQCD, is 
enhanced by the mixing of the J/ip with a vector glueball O that decays preferentially to pn and other VP channels. 
The vector glueball O is required to be fairly narrow and nearly degenerate with the J/t/j. The BES has searched for 
this hypothetical particle in a pn scan across the J/^ region in e~^e~ annihilations as well as in decays ip{2S) -^ ttttO, 
O -^ pn, and found no evidence for its existence |27y2S|] . The data constrains the mass and width of the O to the 
range | mo —'mj/^ \< SOMeV and 4MeV < Tq < 50MeV [Q. This mass, as indicated in Ref. |2^, is several hundred 
MeV lower than the lightest vector glueball observed in lattice simulations of QCD without dynamical quarks. More 
recently, a few more experimental facts unfavorable to this model have been reported by BES. One is the identification 
of isospin- violating VP mode '0(25') -^ wvr^ with a large branching fraction |^. This contradicts the essence of the 
model that the pattern of suppression is dependent on the spin-parity of the final state mesons. The other is the 
finding of suppression of ip{2S) decays into vector plus tensor (VT) final states ||I^. Since hadronic VT decays, 
unlike the VP decays, conserve HHC, some other mechanism must be responsible for this suppression in the model. 
Furthermore, it has been argued that the O may also explain why J /ip decays to (/)/o (named previously S*) but not 
to p5, since the O mixes with the </) and enhances a mode that would be otherwise suppressed [ p6[ . However, the 
observation of non-suppressed ip{2S) — > (pfo Q, which implies the absence of anomalous enhancement in J/ip — > (j)fo, 
would rule out such an explanation. Anselmino et al. extended the idea of J/V' — O mixing to the case of r]c — >VV 
and pp pO| . They suggested that the enhancement of these decays can be attributed to the presence of a trigluonium 
pseudoscalar state with a mass not far from the rjc mass. So far no experimental data has supported the existence of 
such a state. 

Recently Brodsky and Karliner proposed the existence of intrinsic charm | qqcc > Fock components of the light 
vector mesons as another mechanism to account for the J/^ decays to VP channels and their suppression oitp{2S) pOl . 
They also suggested comparing branching fractions for the rjc and ric(2S) as clues to the importance of rjc intrinsic 
charm excitations in the wavefunctions of light hadrons. However, the BES observation of ip{2S) — > unr'^ would again 
appear to disfavor this model. 

Chaichian and Tornqvist suggested a model which invokes a form factor falling exponentially with the energy to 
suppress all "0(25*) decays to lowest-lying two-body meson final states [Q. However, the BES report on observation of 
a number of 0(25') hadronic two-body decays such as bin, cftfo, Ki{1270)K, and ujtt'^ have proved the contrary [|l8|,p[. 
In addition, the BES upper limit at 90% C.L. 5(0^(25) -^ pn) < 2.8 x 10"^ |^ is well below the branching fraction 
predicted by this model, 7 x 10~^. 

The generalized hindered Ml transition mechanism proposed by Pinsky |g^ relates the process 0(25') — > 777' to 
the hindered Ml transition 0(25') — > jrjc- This predicts Q-yri' to be 2 x 10"'^, which, as already shown in Ref. p9[ , 
falls more than an order of magnitude below the BES data (3.6 ± 0.9) x 10~^. According to this model, the hadronic 
decays of 0(25) to VP final states are also a generalized hindered Ml transition. The branching fraction for the decay 
of 0(25') -^ pn is estimated to be 4 x lO"'^, as compared to the measured limit of 2.8 x 10~^. Moreover, it is inferred 
from this model that 0(25) — > 7/2 decay should be suppressed whereas 0(25) — *■ u;f2 should not [g^. However, the 
experimental facts from BES contradict this assumption [pTp] . 

Karl and Roberts have suggested a proposal to explain the pn puzzle based on the mechanism of sequential quark 
pair creation p3| . Even though their predictions could generally accommodate the data for decays of J/0 and 0(25') 
to pn or to K*K, it seems hard to explain the large branching fraction for decays to pn 0] due to the fact that 
their fragmentation probability tends to zero as the mass of the pn decaying system approaches IGeV. 

More recently, Li et al. pi] pointed out that final-state interactions in J/0 and 0(25) decays give rise to effects 
which are of the same order as the tree level amplitudes, and may be a possible explanation for all the observed 
suppressed modes of 0(25') decays including pn, K*K and w/2. They thus predicted qualitatively large production 
rates of aip and K^K for 0(25*), the verification of which may give further support to their model. So far, BES 
has never reported such measurements; nevertheless, useful information on aip and K^K could be obtained from its 
published data as shown in Ref. p7| . The lack of evidence within the invariant mass distribution plots (see Fig. 3 and 
Fig. 5 of Ref. p7J) that the pn recoiled against a p for events of 0(25) — » p^p'^n^ and that n^K^ recoiled against a 
K*'^ for events of 0'(25) — > n^n^ K^K~ suggests that they are unlikely to be the favored modes in 0(25') decays. 

A model put forward by Gerald and Weyers entertains the assumption that the three-gluon annihilation amplitude 
and the QED amplitude add incoherently in all channels in J/0 decays into light hadrons, while in the case of 0(25') 



decays the dominant QCD annihilation amphtude is not into three gluons, but, via a two step process, into a specific 
configuration of five gluons fS^]. Besides explaining the measurements on ip{2S) decays to pi:, K*K, and W7r°, this 
model predicts a sizeable ip{2S) -^ (tt+tt" or r/)/ii(1170) branching fraction. Indeed the BES has performed extensive 
analysis of decays ip(2S) -^ TTTrpn to look for new particles; however, it is unlikely that a conclusive signal for /ii(1170) 
has ever been observed in the inclusive spectrum of ip{2S) decays to ttttptt ||2q| . 

Chen and Braaten proposed an explanation |29|] for the pir puzzle, arguing that the decay J /"ip — > pn is dominated 
by a Fock state in which the cc pair is in a color-octet ^5*1 state which decays via cc — > qq, while the suppression of this 
decay mode for the il'{2S) is attributed to a dynamical effect due to the small energy gap between the mass of the ip{2S) 
and the DD threshold. Using the BES data on the branching fractions into pir and K*K as input, they predicted the 
branching fractions for many other VP decay modes of the ip{2S). Most recently, Fcldmann and KroU parametrized 
the strong interaction mechanism for the hadron-helicity non-conserving decays in a similar way, but interpreted it 
differently |2^. They argued that, for these processes, the charmonium state decays through a light-quark Fock 
component by a soft mechanism, which is characteristic of OZI-rule allowed strong decays. Estimating the light-quark 
admixture by meson-mixing, they also obtained a numerical description of the branching fractions for many VP decay 
modes of the J/ip and ip{2S). The predictions of both models are in good agreement with the measured branching 
fractions (some are preliminary) from the BES experiment Q as well as the PDG data Q. Chen and Braaten's 
proposal also has implications for the angular distributions for two-body decay modes of ■0(25'). Nevertheless, such 
measurements would be extremely difficult, if not impossible, to perform for those strongly suppressed decay modes 
in ip{2S) decays. Feldmann and KroU, on the other hand, has extended their mixing approach to the r]c — > VV decays 
and obtained a reasonable description of the branching fractions for these decays while the 77c (25) -^ VV decays are 
expected to be strongly suppressed. 

From the above discussion we see that essentially none of the models are able to explain all known experimental 
results; in particular no analysis on the suppression ofip{2S) — > VT decays has been given. Not a few models appear to 
have more assumptions than predictions, not to mention quantitative predictions. While the current data seem to rule 
out convincingly some of the models, a few other models may warrant further consideration; for them both detailed 
theoretical analyses and additional experimental tests are demanded. It seems to us that a key premise for physical 
considerations is to establish whether the J/4' decays or the ip{2S) decays are anomalous. An amplitude analysis 
made for the two-body decays of J/-0 to VP [g4| has shown that nothing anomalous is found in the magnitudes of 
the three-gluon and one-photon decay amplitudes. If this is sustained, those arguments presupposing the J/tp as the 
origin of the anomaly should be disregarded. 

In summary, we have examined the pi: puzzle of J/ip and tp{2S) decays in the light of current experimental data. 
The estimates of the ratio of ijj{2S) to J/ip hadronic decay rates, using two different approaches, differ substantially 
from each other. The one using only the data of exclusive hadronic decays appears to be underestimated. The other 
estimate, that is based on computation of the inclusive ggg decay rate by subtracting other decay rates from the 
total decay rate, differs by 2cr-error of this estimate from the naive prediction of PQCD, even though its central 
value is about a factor of two as large as the latter. By comparing this estimated ratio, and taking account the 
corrections associated with specific exclusive decay modes with the corresponding experimental data, anomalies in 
J/ip and 0(25) decays to VP, VT, AP and some other final states are evident. We found from our analysis that the 
exclusive hadronic decays of the 0(25) so far reconstructed experimentally account for only a small fraction of the 
total 0(25) decays and a ratio of ip{2S) to J/-0 hadronic decay rates that exceeds the estimated value is expect to 
occur for a considerable part of the remaining 0(25) decay channels. We have also shown that the recent new results 
from the BES experiment provide crucial tests of various theoretical models proposed to understand the pi: puzzle. 
Further experimental and theoretical efforts are required, in order to fill the missing data and definitively solve the 
perplexing pn puzzle. 
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Table 1. Experimental data on branching fractions for electromagnetic decays of J/'0 and tp{2S) used in our analysis. 
All data are taken from PDG except the branching fraction for ^p{2S) -^ t^t^ which is a first measurement by 



BES |11|. 



Channel 


B(J/V) 


B(^(25)) 


7* — > hadrons 


(17.0 ±2.0)% 


(2.9 ±0.4)% 


e+e~ 


(5.93 ±0.10)% 


(8.8 ±1.3) X 10-^ 


H+fi- 


(5.88 ±0.10)% 


(1.03 ±0.35)% 


T+T- 




(2.71 ±0.70) X 10-^ 



Table 2. Experimental data on branching fractions for J /ij) and i^{2S) decays to lower mass charmonium states 
used in our analysis. All data are taken from PDG B. 



Channel 


B(J/^) 


B(^(25)) 


IVc 


(1.3 ±0.4)% 


(0.28 ± 0.06)% 


w+n-J/tp 




(31.0 ±2.8)% 


tt'W/V 




(18.2 ±2.3)% 


ryJ/V- 




(2.7 ±0.4)% 


ttV/V' 




(9.7 ±2.1) X 10-'' 


7XcO 




(9.3 ±0.9)% 


IXcl 




(8.7 ±0.8)% 


lXc2 




(7.8 ±0.8)% 



Table 3. Branching fractions for the J/ip and ipi^S) exclusive hadronic decays used in our analysis. 
from PDG [1. 



All data are 



Mode 


B(J/V) 


B(V(25)) 


tt+^-tt" 


(1.50 ±0.20)% 


(8 ±5) X 10-^ 


2(^+7r-)7r" 


(3.37 ±0.26)% 


(3.0 ±0.8) X 10-^ 


3(^+7r")7r" 


(2.9 ±0.6)% 


(3.5 ±1.6) X 10-^ 


K+K-TT+Ti- 


(7.2 ±2.3) X 10-=^ 


(1.6 ±0.4) X 10-'^ 


ppn^'K^ 


(6.0 ±0.5) X 10-^ 


(8.0 ±2.0) X 10-4 


ppK^ 


(1.09 ±0.09) X 10-^ 


(1.4 ±0.5) X lO-'i 


PP 


(2.12 ±0.10) X 10-^ 


(1.9 ±0.5) X 10-4 


K+K- 


(2.37 ±0.31) X 10-4 


(1.0 ±0.7) X 10-4 



Table 4. The ratio of branching fractions of ip{2S) and J/ip exchisive decays: Qf — B{il){2S))/B{J jij]). All data 
are from PDG [0. Upper limits are given at the 90% confidence level. 



mode 


HHC 


orbital momentum 


pred. Qf (%) 


meas. Qf (%) 


PP 




1 


24.0 ±5.6 


9.0 ±2.4 


p-K 


0.71 


1 


17.0 ±4.0 


<0.65 


K+K''{892)- 


0.71 


1 


17.0 ±4.0 


<1.1 


c^/2(1270) 




0.71 


17.0 ±4.0 


<4.0 


pa2(1320) 




0.71 


17.0 ±4.0 


<2.1 


/s:*(892)0i?2(1430)" 




0.71 


17.0 ±4.0 


< 1.8 


0/^(1525) 




0.71 


17.0 ±4.0 


<5.6 


6f7r=F 




0.71 


17.0 ±4.0 


17.3 ±5.2 


K:^{1270)K^ 




0.71 


17.0 ±4.0 


>33.3 


K^{UOO)K^ 




0.71 


17.0 ±4.0 


< 8.2 


7?7'(958) 




1 


24.0 ±5.6 


3.5 ±1.0 


IV 




1 


24.0 ±5.6 


< 10.5 



